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Abstract 

We develop a model based on the fractional exclusion statistics, applicable to 
disordered interacting particle systems. Here the species represent elemen- 
tary volumes in an (s + l)-dimensional space, formed by the direct product 
between the s-dimensional space of positions and the quasiparticle energy 
axis. The model is particularly suitable for systems with localized states. 
Using the fermionic, bosonic and Wu perspectives of the FES formalism, we 
present the properties of systems with repulsive screened Coulomb interac- 
tions. We analyze the case of the homogeneous system, as well as a few 
test cases with non-uniformities in the local density of states. The particle 
density is determined and the margin (charging) effects are pointed out. In 
addition, peculiar deviations observed in the temperature dependence of the 
heat capacity and entropy are found, in accordance with the several degrees 
of disorder considered. 



1. Introduction 

The concept of fractional exclusion statistics (FES) was introduced by 
Haldane [H and the statistical properties of FES systems were calculated 
independently by Wu ^ and Isakov jsj. 

The FES was applied to quasiparticle excitations at lowest Landau level 
in the fractional quantum Hall effect, spinon excitations in a spin-| quantum 
antiferomagnet |2|, Bose and Fermi systems described in the thermody- 
namic Bethe ansatz [sl, 0, [s], 0] , excitations [H or motifs of spins 0, Isf in spin 
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chains, elementary volumes obtained by coarse-graining in the phase-space 



of a system g, i, El, [III [Til, [l3 



etc. 



More recently, a stochastic method for the simulation of the time evo- 

as a generalization of a similar 



lution of FES systems was introduced [14 



method used for Bose and Fermi systems [15|, and the general FES formal- 
ism was amended to include the change of the FES parameters at the change 
of the particle species ll|, |l6|, ll7|, [iSj . This makes of FES a general formal- 



ism which is able to describe interacting particle systems as ideal gases of 



quasiparticles [19|, l20 



A FES system consists of a countable set of species, indexed by i,j = 
0, 1, 2, . . .. Each species contains a finite number of single-particle states and 
particles, denoted by Gi and Ni, respectively. The number of states in the 
species depend on the number of particles. For small variations of the number 
of particles around some reference distribution, {A'j}j=o,i,...; the number of 
states changes by 

= -^a./AT^-, (1) 

j 

where by (5iVi we denote the particle variations and Oj/s are called the FES 
parameters 



The FES parameters must satisfy certain rules [16|, namely if we split an 
arbitrary species, j, into a number of sub-species, jo,ji, . . ., then all the pa- 
rameters aki, with both, k and / different from j, remain unchanged, whereas 
the rest of the parameters must satisfy the relations: 



a 



a 



aijo = ttij^ = . . . , for any i, j 

aj^i + aj^i + for any z, j 



These rules are satisfied by the ansatz [13 



(2a) 
(2b) 
(2c) 



(3a) 



where the parameters Cifj-, called the "extensive" parameters, are propor- 
tional to Gi 

(3b) 



O-ij — OjijKjxi. 



The parameters af^ refer to only one species and do not depend on Gi. 
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The number of microscopic configurations compatible to a given distri- 
bution of particles on species, {Ni}, is 



m{G,-i)\ 



and W^{{N,i}) 



m{G,-Ni)i 



(4) 



if the particles are bosons and fermions, respectively. 

If for each species of particles, say species i, we associate an energy, 
Si, and a chemical potential /ij then the equilibrium particle distribution, 
rij = Ni/Gi, is obtained by maximizing the partition function. 



Z(^) = ^ jiy(^)({iV,})exp 



(5) 



with respect to the distribution {Ni}, taking into account that the G/s vary 
with {Ni} according to ([1]). In Eq. (|5]) P is used instead of "B" or "F" from 
Eqs. (jl]) and f3 = 1/(/cb7') is the inverse temperature. 

The maximization of Z^^^ with the conditions ([T]) gives 



1 ih 

- e,) + In ' = ± Z-/ ^^[^ ^ 

j 



Hi 



(6a) 



where the all the upper "+" signs are for bosons and the lower "— " signs are 
for fermions. 



If one uses the ansatz ([3]), Eq. (l6al) becomes [13 



In 



PM-a 



± Gjttji ln[l ± n^]. 



(6b) 



In FES, a system of fermions with a set of parameters, {oij}, may be 
interpreted as a system of bosons with the parameters {oij + 6ij} and vice- 
versa, a system of bosons of parameters {oij} may be interpreted as a system 
of fermions with the parameters {aij — 6ij}. Therefore it is more natural to 
refer to the Bose and Fermi formulations, rather than to bosons and fermions. 

The most used formulation of FES is the one employed by Wu in Ref. jij, 
which will be denoted here by "W". To see how this is related to the Bose 
formulation we define the number of states in the absence of particles in the 
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system, = Gi + Yl^ij-^j^ ^"^^ ^ particle population, = Ni/Gf^. 
Then the nj^'s are determined in two steps. First one solves the system 

to determine the Wj's, and then the ra^'s are calculated from 

J2i^^,w, + a.,GY/GT)n, = 1. (7b) 



Comparing Eqs. f[6al) and fITal) we observe that Wi = 1/nf 



In Refs. [19|, |20[ one of us have shown that FES is the natural description 
of an interacting particle system as a system of ideal quasiparticles. In this 
paper we exemplify this procedure on disordered systems of interacting lo- 
calized particles. This method may find applications in glassy systems such 
as the Coulomb glasses and systems of bosons trapped in optical lattices. 

The structure of the paper is as follows. In Section H] we introduce our 
model, in which species are elementary volumes in the (s + 1) -dimensional 
space formed by the direct product between the s-dimensional space of po- 
sitions, Q, and the energy axis, e or e. These species are related by FES 
parameters, which we calculate. The formalism is then applied to homo- 
geneous systems (mean-field approximation) and the FES parameters are 
calculated for a repulsive screened Coulomb potential. In Section 12] we ap- 
ply our formalism on a few test cases, which are one-dimensional periodic 
systems with several degrees of disorder. 



2. Model and formalism 

The particles are localized on random sites in a solid s-dimensional ma- 
trix. The positions of the sites are denoted by r/, J = 1, 2, . . . , A^o, where 
No is the total number of sites. We assume that the wavefunctions of the 
particles do not overlap and the total energy of the system is 

E = ^^vjnr, + IY1 ^rirj^ririrj . (8) 

where is the energy and rij-j is the occupation number of the site /; the 
total particle number is = n,.^. 
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We shall work in the continuous limit, so we define the density of sites, 
cr(r, e) = XI/ '^''('^ ~ ^/)^(^ ~ ^r/), and the particle density, p(r, e) = 
ri)6{e—erj)nrj. The average particle population in the (s+l)D volume 6Qx6e 

is defined as n(r, e) = [Js^d^^ Ise^^ P^^^^^yi.Isn^'^^ /se '^^ '^(''' ^)]' where we 
assume that the volume is large enough, so that J^^ d'^r J^^ dea{r, e) > 0. In 
these notations the total energy of the system ([HD becomes 

E = d'r ep(r, e) de + - / d'r / d'r' 

/*€max /'^max 

X / de rfe'p(r,e)p(r',e')Kr' (9) 

^min *^ ^min 

where Vt is total the volume of the system and [emin, ^max] is the interval in 
which e takes values, with emin > 0. 

We shall assume that emin = 0, e^ax ^ ksT (where /c^ is the Boltzmann 
constant and T is the temperature) and the interaction energy depends only 
on the distance between the sites, i.e. K/rj = ^(k/ ~ Because emax S> 
ksT, we shall take emax = oo. 

For concreteness, we analyze only Fermi system, but the formalism can 
be easily extended to bosons. 

To apply FES, we have to divide the system into species. We do this by 
coarse-graining the parameters space of the system, f2 x e, into the elementary 
volumes, 5VL^ x 5ei. By the lower case Greek letters, e.g. ^ = 0, 1, . . ., we 
identify the spatial volume and by the lower case Latin letters, e.g. i = 
0, 1, . . ., we identify the energy intervals, 5ei = [e^, Cj+i]. We have eo = 0. 

We identify a species either directly, by 5Vt^ x 5ej, or by the subscripts, 
We assume that each elementary volume, 5Vt^, is centered at and 
contains a large enough number of sites and particles to justify the appli- 
cation of the statistical methods and, in particular the Stirling approxima- 
tion for the logarithms of factorial numbers. In each of the volumes, say 
we have a distribution of sites, r/ G 5VL^, of energies e^j. Under the 
assumption that the number of sites is large enough, the set of energies 
{^r/jr/e^Vj form a (quasi) continuous distribution along the e axis, with a 
density (J^{e) = J^^^ a{r,e)d^T. The number of states in the species (^, z) is 
then G^i = a^{ei)6ei and the number of particles is N^i = a^{ei)6ein{T^, e^) = 
SQ^Seip{r^,ei). We shall use the notations p^(e) = p(r^,e) for the particle 
density and = f dep^r^, e) for the volume particle density. 



We define the quasiparticle energies, er^, in a similar way as in Ref. [12 
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l2li,by 



(10) 



Because of the identity E = ^nfirj, the thermodynamics of the quasipar- 
ticle gas follows identically the thermodynamics of the original system. 
In the continuous limit, Eq. (fTO|) becomes 



d^'r / (ieV(|r/ - r|)a(r, e)?7.(r, e) 



(11) 



which is similar to Eq. (5) of Ref. 12|, except that here we have an extra 
integral over r. 

We have a new parameters space, f2 x e, in which, by construction, e^i^i = 
emin = and e^ax = ^max = oo- Using Eq. (fTTj) we obtain the new DOS, 



a[r,e(e)] = a(r, e) 



di 
Te 



-1 



o- r, e 



|l + J^d'r'V{\r - r'|)a(r',e)n(r',e)| ' 



(12) 



Assuming that Eq. (fTTl) defines a one-to-one function e(e) - which may 
also be inverted to e(e) - we split the space fixe into species, x 6ei, as 
we did with fi x e, in such a way that ii = e(ej) for any i and each species 
contains G^t states and N^t particles, as before. 



By direct application of the procedure of Ref. [12, l21| to this interaction 
potential and particle species we obtain the FES parameters 



SijC^ii^i) + 9{i-j)6e 



dei 



6ij6Q^a{ri:, e^) + 6{i - j)6ei6Q. 



V{h-r,\) 



Vih-r,\) (13) 



where the first doublet, (^i), specifies the species in which the number of 
states changes, whereas the second doublet, (r/j), specifies the species in 
which the number of particle changes; 6{k) is the step function, 6{k > 0) = 1 
and 6{k < 0) = 0. The manifestation of the FES parameters given by Eq. 
f|T3|) for a system with cr(r, e) = cr(r), i.e. independent of e, between species 
of the same quasiparticle energies is represented in Fig. [H 

We observe that the FES parameters ( fT3l) obey the rules ([2]) IGj and 
define a new ansatz, which is a generalization of the ansatz ([3]) jl3| . 
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2.1. Equilibrium thermodynamics 

Since we have fermions in the systems, we employ the Fermi formulation 



13[ to calculate the equilibrium thermodynamics. Plugging the a parameters 



f lT^ into the equations we get 

= - 6,) + In ^ 'J"/^""!''^ ^ + ^ a,,;e. ln[l - n^(r„ g,)]. (14) 

In the continuous limit, Eq. ( |T^ becomes 

= /3(^-e) + lnl^:^^^^+ /"rfVV(|r'-r|)(T[r',e(e)]ln[l-n^(r',e)] 

ln[l-n^(r',e')]. (15) 



n 



oo 



+ / dVy(|r'-r|) / d~e' ^""^""'''^ 



9e 



B 



In the bosonic formulation, as we mentioned above, a^i-^j = S^rt^ij + a^i-r^j, 
and the dimension of the species is the number of available states, = 
G^i — N^i + 1 (we used the notation G^^ to avoid confusion with the number of 
states in the species, G^i). This changes also the definition of the population 
to nP{r^, ei) = N^i/G^^ = n{r^, ej)/[l — n(rg, ii)]. Plugging the new quantities 
into Eqs. we obtain the system of equations for n^{r^, ei): 



= /3(/.-6.)+lni±^^i^-5^<,,ln[l+n^(r,,e,)] (16) 



m 

and in the continuous limit 



"1 



= ^(^-f) + ,ni±!Li£i£) 



- I rfVy(|r'-r|) / d~e' 



f d'r'V{\r' - r|)a[r', e(e)] ln[l + n^(r', e)] 
Jn 

ln[l + nB(r',e')]. (17) 



de 



e(e') 

In Wu's formulation Q the dimension of the species is G^ = G^i + 
Er,i oiiKvj^vj = + ^rij Oil;vj^vj ^nd n^{r^, e^) = N^i/G"^. The equilib- 
rium population is calculated from 

(1 + w^,) n (j^) = e(^^>-'^)/'=^, (18a) 

-L-L \1+Wrijj 



and 



{Siir,SijWr,j + (3^i,r,j) U^j =1, 



(18b) 



where P^i^nj = = ^(v^ij + '^(iivj^vj ■ observe again that 

Having the quasiparticle populations, we can calculate any thermody- 
namical quantity. The internal energy of the system is 

/ £,i e.i 

(19) 

or, in the continuous limit. 



U{T, /i) 



d'r / rfep(r,e)e, (20) 



where P = F, B, W and 5"^(r, e) is given by Eq. (fT2|) . The other two 
densities of states are 5"^(r, e) = 5-^(r, e) — p(r, e) and ^^(r, e) and all satisfy 
the relation p(r, e) = cr^(r, e)n^{r, e) = (T^(r, e)n^(r, e) = 5"^(r, e)n^(r, e). 

We shall use the notation p^{e) = p(r^,e) for the quasi-particle density 
in fl^, at quasiparticle energy e, and = J dep{r, e) for the volume quasi- 
particle density. One should note that p^ = p^ by definition. 

In order to define the DOS ^^(r, e) we have to define the densities of the 
FES parameters [Ul, ll3|, ttre-r'e', by 



where 



(21a) 



mr5-r,|)(21b) 



Using Eqs. fl2Tl) we write cr^(r, e) = a^(r, e) + d'^r' de' ari-yi'Pi^', e')- 
Similarly, the total particle number is 



N{T, p) 



/ d'r / dea^{r,i)n^{r,i) = / d'r / dep{r,i). (22) 
Jq Jo Jn Jo 



The heat capacity and the entropy of the system are 



[df 



dU{T, p) dU{T, p) dN{T, p) ( dU{T, p) 



N 



dT 



dp 



dT \ dp 



(23) 



(a) 



e. 



(b) 



'5 



Figure 1: Partitioning of real space and quasiparticle energy axis. Upon inserting an extra 
particle in species i), the numbers of available states in all species (77, i) are changed. 



and 

S = kB \n{W^ 
respectively, and they satisfy the equation 

'dS_ 



T 



(24) 
(25) 



N 



2.2. Homogeneous system 

If cr(r, e) = o-Q{e) is independent of r and if we impose periodic boundary 
conditions or ignore the effect of the surfaces (deep inside the sohd), then 
both, n(r, e) and 5"(r, e), are independent of r and in the Eqs. (fTSl) and f|T71) 
we can perform the integrals over Q to obtain an effective FES parameter. 



aefr(e) = an[e{e)]h 



de 



de 



(26) 
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where 

Iv= [ rfVr(|r'-r|). (27) 
Jn 

Plugging Oeff back into Eqs. (fTSl) and (fTTI) we get the equations 

- ^) ^ - m I' - -;<■-■ - "'b';- "r°°" (28) 

ra^(r,e) n^(r, e) 

3. Applications 

We exemphfy the formahsm above by assuming a two-particle screend 
potential of the general form 

V(r;7.A) = K5^H<^. (29) 

where r = |r— r'|. In particular, if 7 = 1, we have the usual screened Coulomb 
and Yukawa type potentials. Such systems exhibit standard thermostatistical 



behavior if (i/7 < 1 [22| and the interactions are classified as short ranged. 
If (i/7 > 1 the systems obey non-extensive thermodynamics, i.e. quantities 
like total energy are not extensive due to the long range interactions. 

To remove the singularity at the origin that appear in the integrals over 
V in Eqs. (fTSl) . (fT7|) and (1271) . we introduce a cut-off at radius Rq, below 
which the potential remains constant - V{r) = V{Ro) = Vq for r < Rq. With 
this assumption, for a homogeneous system. 



27rt 



r (!) Jo 



POO 

/ dr r^"V(r) = lyi + Iv2, (30) 
Jo 



where 



Ivi = / dr r^- V(r) = —. ^R'Vo 

T{^)Jo r(f + i) ° 



/y2 = Iv-Ivi = ^^^y-'r(^s-^,^y (31) 
Here r(s,a;) = t^~^e~*dt is the upper incomplete gamma function. 
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For a screened Coulomb-type interaction in a s-dimensional system (s 
1, 2, 3), the term Iy2 can be expressed as: 



4^ = 2KaEi(i?o/A) 

t2D o \ f 

Iy2 = ZTTKcrAexpl — — 
4f = 47iKa \{Ro + A) exp , (32) 

where Ei{z) is the exponential integral, 

poo -t roo -zt 

E^{z)= / — dt= / dt. (33) 



5". i. Application to one- dimensional disordered interacting systems 

To illustrate some peculiar thermodynamical properties of interacting dis- 
ordered systems we apply our FES formalism to ID disordered systems with 
periodic boundary conditions. The physical space is partitioned into Nr sub- 
volumes on the coordinate axis and Ne segments on the energy axis, which 
define a total number of Nn x Ne species. The periodic boundary conditions 
are imposed using the minimum image convention, i.e. the interaction drops 
at a distance equal to half the length of the repetitive unit. 

For simplicity we shall assume in the following that (T(r, e) [and cr^(e)] are 
independent of e, so we shall simplify the notation of the DOS to cr(r) (and 
0"^). In this case the FES parameters f[T^ simplify to 

o^ii-m = %crc(ei)^(|r5 " ^r,!) = ^^^^^^^^(r^, ei)V{\r^ - (34) 

which is diagonal in the indices i and j. 

For our ID system we apply a spatial partitioning such that = = 1 
and we take the cut-off distance, Ro, to be half of the species dimension, 
i.e. Rq = -R5/2. Furthermore, we set V{R^) = 1, which implies V{Ro) = 2. 
Under these assumptions we rewrite the FES parameters in Eq. 0341) as: 

= o-^{(^i)VoR^ 
o;^i;nj = Sij(T^ie-i)iVo/2)Rl/\r^-rn\ for i ^ j. 

The quasiparticle densities p^(e), can be found in any of the descrip- 
tions aforementioned [F, B, W) by solving the corresponding nonlinear sys- 
tem 0151) . f[T7|) or (ITSl) . The solution is found iteratively using gradient descent 
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- as implemented by GSL - starting from an initial guess solution. Since the 
thermodynamical quantities of interest should be obtained for a given par- 
ticle number rather than for a given /i, an extra loop is introduced so 
that the corresponding chemical potential is found in agreement with the 
normalization relation fl22|l . 

We consider a system with repulsive screened Coulomb interactions, pa- 
rameterized by 7 = 1 and A = 3R^. For the partition of the physical space 
and the coarse graining on the energy axis we take Nr = 20 and A'^^; = 50, 
respectively. The system contains a number of N = NuEpa particles, where 
the Ep = 1 is the Fermi energy. In this way we set up the energy scale. We 
also assume that /c^ = 1. 

Two particularly transparent examples of an inhomogeneous system are 
obtained by introducing a dip or a peak in the density of states of a ho- 
mogeneous system. More concretely, we pick one sub-volume C, for which 
the density of states is half or twice the reference (constant) value of the 
other sub- volumes, respectively: = 1.0 for any C, ( and cr^ = 0.5 ("dip" 
case) or 2.0 ("peak" case). For symmetry reasons we have considered Nji + 1 
sub- volumes. Figure [2] depicts the particle densities for the two types of sys- 
tems considered at T = 1. Unlike a homogeneous system, where there is 
a unique function n^i = rii for all sub-volumes, corresponding to a (single) 
mean field FES parameter agg, we now have different populations with a 
spatial distribution. As one can see from Fig. [2] the systems obey a mirrored 
symmetry: the density of particles in sub-volume C, drops in the "dip" case 
and is enhanced in "peak" case, as compared with a mean field calculation for 
a constant a = 1.0. The particle densities in nearest neighbors sub-volumes 
of (, i.e. C ~ 1 ^^'^ C + 1) exhibit deviations from the mean field values and 
opposite to the values in sub- volume (. This can be explained as follows. 

At equilibrium, in the case of a dip in the density of states, the repulsive 
interaction between particles drives accumulations of particles towards the 
edges of the region with constant a = 1.0, as the number of particle in sub- 
volume ( is lower due to a smaller number of available states. The opposite is 
found in the case of a peak in the DOS, where the larger number of particles 
in sub- volume ( repel the particles in the adjacent regions. Starting with the 
C — 2, C + 2 sub- volumes, the populations already get very close to the mean 
field values. 

However, for dips or peaks of finite widths, a broader distribution of 
particle densities can be observed in the insets of Fig. [21 Here we took the 
linear dependences 0"^+^ = = 0.5 -|- O.lr (dip) and 0"^+^ = cr^_^ = 
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2.0 — 0.2r (peak), with r = 1, 2, 3, 4. Note that due to the symmetry, except 
the ( sub- volume, we have the pairs P(-r,i = Pc+T,i, for any r > 0. 

We further analyze the particle distribution in a periodic one-dimensional 
system with an interface between high [a^ = 1.0, for ^ < Nji/2) and low 
(a^ = 0.5, for ^ > Npi/2) DOS regions. The obtained results are represented 
in Fig. [3] for a temperature T = 1. From the main plot one can see the 
particle densities are divided in two groups corresponding to the two values 
of cr. Like in the previous case, significant deviations occur in the sub- volumes 
which form the interfaces, which can be understood by a similar reasoning. 
The particle densities in the other sub-volumes are depicted in detail in the 
two insets, indicating a convergence towards the values corresponding to the 
mid-points of the two constant a regions. 

Now we analyze the properties of a system with disorder which we intro- 
duce by randomly distributed values for the local density of states at each 
site. More specifically, we consider a step distribution, centered around the 
average value cxo = 2 and several degrees of disorder are set by the width of 
the distribution step Acr. Figure IH^a) shows a typical particle distribution 
obtained at a temperature T = 1 and maximum disorder. Act = 4. 

Like in the two previously analyzed examples one can observe the parti- 
cle densities are larger at sites with larger DOS, although now we find the 
disorder specific distributions p^j. In the lower plot is represented the local 
DOS together with the particle distributions at each site ^, for two temper- 
atures Ti = 0.5 and = 6. One can see that for the higher temperature 
Th the particle distribution follows closely the distribution, while for the 
the lower temperature T/ the particles are more evenly distributed. In the 
high temperature limit the interactions become less and less important, the 
particle distributions approach the Maxwell-Boltzmann distribution and all 
available states become equally probable. Consequently, the spatial particle 
density in each sub- volume becomes proportional with the density of states 
a^. By contrast, in the lower temperature limit, the repulsive interactions 
tend to level the particle distribution. At low and intermediate temperatures, 
there is no clear relation to the local DOS anymore and the exact equilibrium 
distribution can be obtained using the FES formalism. 

The results presented so far for the particle density distributions may 
be obtained in any of the three descriptions mentioned in the section [2l i.e. 
the fermionic, bosonic and Wu pictures, which all lead to the same physical 
results. To further illustrate this fact, we plotted in Fig. |5]the populations 
and the quasi-particle density of states in the fermionic and bosonic pictures. 
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for the same disordered system analyzed before (ctq = 2, Act = 4, T = 1). As 
expected all the populations in fermionic description, nj-, are smaller than 
unity, while some of the bosonic populations, n^^, indicate a large increase at 
low energies. This difference regarding the populations in the two pictures is 
compensated by the behavior of the quasi-particle density of states, and 
a^, so that the particle densities p^i and are the same in any description. 
Specifically, we have for low energies < a^, while at high energies both 

and asymptotically approach the Wu density of states, . 

The heat capacity of such disordered systems exhibits peculiar effects. In 
a uniform system, i.e. with a constant DOS, the heat capacity is the same for 
any diagonal FES parameter a 23,^. However, here we have FES parame- 
ters of the form (IT^ . i.e. diagonal in energy and non-diagonal with respect to 
position indices. The results are presented in Fig. M^a) for different distribu- 
tions of the local DOS, as indicated in the inset, corresponding to the values 
Act = 1, 0.75, 0.5(Jo. By increasing the disorder a deviation from the reference 
constant-DOS dependence of the heat capacity per particle is observed, with 
a maximum above 1. Similar deviations are observed in the temperature 
dendence of the entropy in Fig. EJ^b), which are in agreement with relation 
( 125|) . These deviations vanish in the low disorder or high temperature limits. 

We also performed ensemble averages on disorder and a qualitatively 
similar behavior was found. 



4. Conclusions 

We have formulated an approach based on the fractional exclusion statis- 
tics to calculate the thermodynamic properties of disordered interacting sys- 
tems. The interacting particles are localized at randomly distributed sites 
in a solid matrix. By partitioning the real space besides the usual division 
of the energy axis and relating the obtained species by FES parameters we 
are able to go beyond a mean field picture, describe the particularities of 
such disordered systems and calculate their thermodynamic properties in an 
intermediate temperature range. 

We also analyzed systems which do not appear to be disordered at the 
scale at which we partition the system, but exhibit a non-homogeneous dis- 
tribution of sites. We calculate the particle distributions, the quasiparticle 
densities of states and we emphasize the margin effects. 

Beside its fundamental aspect, this method may be applied in the calcu- 
lation of the thermodynamic and transport properties of mesoscopic systems. 
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Figure 2: (Color online) Quasiparticle density, p^(e), for systems with a dip (cr^ = 0.5) 
(a) and a peak {a^ = 2.0) (b) in the otherwise constant density of states (cr^ = 1.0, for 
^ ^ <^). The symbols represent the particle densities which correspond to the species ( 
[squares (a) and circles (b)] and its first nearest neighbors (up/down triangles). The solid 
black lines correspond to a uniform systems with a = 1.0. The insets contain similar data 
for systems with a finite width dip/peak. The dashed lines represent particle densities 
for sub- volumes [C — t, C + t], with r < 4. The volume density of (quasi)particles in 
each species are represented in (c). Pairs of vertical dotted lines mark the extension of the 
considered dips and peaks. 
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Figure 3: (Color online) (a) Quasiparticle density for a system with a step-like density of 
states, i.e. cr^ = 1. for ^ < Nji/2 (upper eurves, in red) and — 0.5, for ^ > A^_r/2 (lower 
eurves, in blue). The symbols (up/down triangles) correspond to the four sub-volumes 
located at the interfaces. In the insets (left, for ^ < Nji/2, and right, for ^ > A^_r/2) are 
shown details of p^(e) for the remaining sub- volumes. The numbers indicate the minimum 
distance from either interface, (b) Spatial (quasi)particle density of the considered 
system. 
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Figure 4: (Color online) (a) Quasiparticle density for a disordered system {(Tq — 2, Act = 
4). The species with high (blue/solid) and low (red/dashed) DOS are indicated, (b) 
Density of states cr^ (black/solid with dots) and particle numbers in each sub- volume for 
a high {Tfi = 6.0) (red/dashed) and a low (T; — 0.5) (blue/solid) temperature. 




Figure 5: (Color online) Populations and quasi-particle density of states (red/dashed for 
cr^ > 2, blue/dashed for < 2) in the fermionic (a-b) and bosonic (c-d) descriptions. 
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Figure 6: (Color online) Heat capacity (a) and entropy (b) for different degrees of disorder: 
Act = (black/solid), 0.5 (red/dotted), 0.75 (green/dashed) and 1 (blue/dashed-dotted), 
as indicated in the inset. 



20 



